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I ■ We discuss the angular momentum blockade in small d-wave superconducting 
p ,! 

^ . grains in an external field. We find that abrupt changes in angular momentum state 
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of the condensate, angular momentum blockade, occur as a result of changes in the 
angular momentum of the condensate in an external magnetic field. The effect rep- 



. resents a direct analogy with the Coulomb blockade. We use the Ginzburg-Landau 



formalism to illustrate how a magnetic field induces a deviation from the d-wave 
symmetry which is described by a [d,^2_y2 + idxy)-OTdei parameter. We derive the 
' behavior of the volume magnetic susceptibility as a function of the magnetic field, 



and corresponding magnetization jumps at critical values of the field that should be 
experimentally observable in superconducting grains. 
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I. INTRODUCTION 



The precise nature of the superconducting state in cuprate superconductors has been 
discussed extensively since the discovery of high-T^ superconductivity jl|. While most of 
the data can be well covered assuming a pure rf-wave symmetry of the order parameter 
2|-[5|, the identification of the precise symmetry of the order parameter remains one of the 
active areas of research. One can imagine that the pairing state symmetry is affected by the 
crystal field and, as is the case of YBCO, by the presence of oxygen chains js, 7|. Moreover, 
even if the pairing state symmetry is a simple d-wave, it can be modified and distorted by 
application of an external field js] and by scattering off defects jsHl. These distortions 



can also depend on the doping and the nature of correlation effects in these materials [12 1. 
One might for example expect that the symmetry of the superconducting order changes 
as a function of doping and therefore this pairing symmetry contains useful information 
about the microscopic interactions responsible for pairing. Indeed recent results suggest 
that nanoscale (i-wave superconductors can be fully gapped and this minimal gap (on the 
scale of 10 mK) can be modified by an external magnetic field 13||. We thus feel that the 
whole subject warrants a fresh look in the light of recent findings. 

We would like to revisit the question of the gap induction by a magnetic field in a 
nanoscale (i-wave superconductor. While the general expectation that a magnetic field will 
induce additional components of the order parameter remains, the specific case of a small 
superconducting grain allows for sharp transitions between states with different orbital mag- 
netic moment carried by pairs. These changes in magnetization can be observable in the 
case of small d-wave grains, as we will point out. Therefore qualitatively new effects can be 
expected in investigating small grains of rf-wave superconductors. Earlier it was pointed out 
by Laughlin that the presence of nodes makes the d-wave superconducting states inherently 
unstable to the induction of novel components of the gap 8|. Similar effects of induction of 
additional components are expected in the presence of a steady supercurrent Q,|l5l. The 
general source of instability of the pure ci-wave gap is the presence of the nodes. A sec- 
ondary component allows these nodes to be completely lifted. The most probable secondary 
components proposed in this context are d + is jl^ and d + id' js] states. Both of these 
states will " seal the node of the gap" . Both of the states break time reversal symmetry and 
thus will induce edge currents. The key difference between these two proposed states is that 
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a d + id' state is chiral and hence will have an orbital moment carried by Cooper pairs. 
The appearance of an id'-wave competing order component is expected in the presence 
externa, magnetic field fl Q. and has been used to exp.a,n some experiments & 
|19I |. In the present paper we will consider values H ~ Hc2 for the external magnetic field, 
such that the high density of produced vortices almost destroys the superconducting state. 
In this region of the parameter space, the action of the external field is able to induce a 
distortion of the original d-wave state into a d^^-y^ + id^y (briefly referred to as d + id') 
bulk state, having an intrinsic angular momentum. That induces a movement across ground 
states with different symmetries, and consequent jumps in thermodynamic quantities. The 
low-energy quasiparticles sitting in the nodes of the (i-wave superconductor (SC) have a 
vanishingly small gap and allow for the generation of a secondary component; this feature of 
(i-wave SCs is called marginal stability, which is prevented in the s-wave SC by the presence 
of a finite gap on the Fermi surface 20|. The changes from the zero angular momentum 
state to a state with finite angular momentum will occur via a set of steps corresponding to 
Lz = 0,2h,4h, .... These steps are small and are not observable in a bulk system since the 
number of Cooper pairs is large and the effect is therefore oc 1/Npairs- Only in a sufficiently 
small grain, this staircase of angular momentum jumps becomes evident. In this paper we 
will analyze the angular momentum jumps in a d-wave grain as one sweeps the applied 
magnetic field. These jumps are similar to the Coulomb blockade charge jumps seen in 
quantum dots 21|, see Table HI The underlying energy that controls these jumps is the 
energy of orbital moments in magnetic fields and thus is much smaller than the Coulomb 
energy controlling charge blockade phenomena. So in analogy with the Coulomb blockade 
we propose that the field-induced angular momentum changes in the Cooper pair states 
represent an angular momentum blockade that can still be seen in the susceptibility for a 
small superconducting grain. This is the main finding of this paper. 



Existing results [22| allow us to exclude a similar blockade effect in ordinary s-wave 
superconductors, in agreement with a vanishing contribution of the angular momentum 
term to their free energy. Indeed, the magnetic susceptibility is therein proven to show a 
maximum for a well-defined value of the magnetic field as soon as ultra small SC grains are 
considered. 

We also point out that the effects discussed in this paper will equally be present in the p- 
wave superconductors. The obvious analogy will be the induction ofp+ip' superconductivity 
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in the presence of a magnetic field due to tlie magnetic moment coupling to an external field. 
We will not specifically elaborate on that case but wish to point this obvious extension of 
the calculations presented. 

The plan of this paper is as follows. First we will connect the field-induced component of 
the order parameter to the angular momentum. This will allow us to express the magnetic 
observables in terms of the angular momentum itself, by adopting the Ginzburg-Landau 
description for the superconductor. We will conclude by analyzing the dependence of these 
observables on the external magnetic field, and by estimating the orders of magnitude in- 
volved in the effect. All the relevant numerical values and some technical details can be 
found in the Appendices. 

II. ANGULAR MOMENTUM BLOCKADE 

First we elaborate on the notion of angular momentum blockade. The low energy gaps 
induced by the magnetic field translate into very long length scales relevant for the formation 
of the new component of condensate. Therefore the formation of the induced gaps can be 
described by a continuum theory. In that limit the total angular momentum of the 
condensate is a good quantum number. Simple inspection shows (see below) that the pure 
c/-wave state has zero angular momentum. The d + id' state is a chiral state and has a finite 
angular momentum ^ 0. Therefore any changes from the ground state with zero quantum 
number to a ground state with nonzero quantum number would proceed as a set of steps in 
Lz- The smallest steps of total angular momentum one can get in the SC grain will be in 
units of 2h. Indeed we envision the set of incremental steps by which Cooper pairs convert 
from dio d + id' as a set of jumps in L^. These jumps result in jumps of magnetization that 
can be seen only if the relative change is large. Hence it is expected that these changes will 
be seen only in small samples. This discussion is analogous to the Coulomb blockade that 
describes single electron charge jumps. Table HI 

To illustrate the mechanism we will use the Ginzburg-Landau theory for an high-T^ SC 
in an external magnetic field, which is written in terms of the complex order parameter 

A(f,e) = exp[z<l>(r)]A(e) , (1) 
where $(r*) is the phase of the order parameter, and 9 refers to the dependence of the 
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pair wave function on the position on the (2D) Fermi surface. The physical origin for the 
generation of the second component id' resides in the induced bulk magnetic moment (M^), 
and the relevant interaction is the {Mz)B coupling to the magnetic field B \\ z, where B is 
the magnetic induction in presence of an external field H. The d-wave state before turning 
on the magnetic field can be regarded as a superposition of the pairs with orbital momentum 
L, = ±2h: 

Ao{e) = Ao cos2e = ^ [exp{2ie) + exp(-2z0)] , (2) 

where Aq is the gap magnitude of the d component, and a 2D geometry of the Fermi surface 
is considered in the model, because of the layered structure of the cuprates. The external 
field H has the effect of shifting these two = ±2h components, linearly and with opposite 
signs as: 

Ao(^) ^ ^ [(1 + vB) exp{2te) + (1 - r]B) exp(-2^^)] = Ao(^) + iB\A^{e)\ , (3) 

where |Ai(^)| oc 77 sin 26' is the modulus of the d' component and 77 is the coupling constant. 
Unlike each of its two components taken alone (having nodes of the gap), the full d + id' 
state is fully gapped [16|, |23| , see Fig. [H 

To a good approximation, the angular momentum of the whole sample about z is 
related to Aq, Ai (see Appendix Rl for derivation) by: 

|Ai| = |Aol /f;^.' . (4) 
2n g {V/a^c) 

Here the components of the order parameter have been assumed to be spatially uniform 
within the sample, which is justified to a good approximation for sample sizes much smaller 
than the characteristic length of variation of the SC wave function; (a^c) denotes the unit 
cell volume of the SC's crystal structure (being c the edge || i || -B field, a ~ 4A, c ^ lOA), 
V the volume of the SC sample, and g ~ 0.1 the superconducting fraction of the material. 

Let us define the non- negative integer variable / = \ {Lz)\/{2H) and the characteristic area 
5*0 = a^c/{d ■ g), where d is the sample thickness and V = Ad is the volume. This gives 



(5) 
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(a) Competing order (b) Full gap 

FIG. 1. (a) The black circle representing the Fermi surface is shown together with the angular 
dependent gap amplitudes. The dashed green curve is the d-wave gap Aq, and the smaller dotted 
orange curve is the induced d'-wave gap Ai. For display purposes, the relative proportion of the 
gaps has not been preserved in the plot, (b) Outside the black circular Fermi surface, the total 
gap Y^|Aop + |Aip of the fully gapped d + state is shown in blue. The former gapless node 
lines of the d-state, which are the diagonals w.r.t. the [100] and [010] crystallographic directions, 
are sealed by the induced imaginary component id! ^ which renders the overall d + id' state fully 
gapped. 



In the case of an isotropic gradient tensor i^jj, and if the grain is small enough to allow 
us neglect the spatial dependence of Ao(r), Ai(r), the GL functional takes the form 



'v=Ad 



a |Ao| 
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(-z-A)Ao 

c 



r/5|Ao||Ai 



iVn 



|Ai| 



(6) 



provided that we omit the magnetic field self-energy term, which does not affect the behavior 
of the magnetization or the susceptibility that we are going to describe. In the Coulomb 



gauge A{r) = ^{-y,x), Eq. 



d^r 



V=Ad 



a lAo 



gives 

2 

4 



c2 4 



Aor + i^^^r^lAol 



r,B |Ao| — / + — |Ao| j^l 

(7) 



The coupling constant rj is derived and discussed in the Appendix of Ref. 161 ]. 
In the same way as the Coulomb blockade leads to jumps in charges and in charging 
energy, the analogous phenomena for the angular momentum blockade involve changes in 



angular momentum and magnetization. There will be jumps and spikes in magnetization 
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COULOMB BLOCKADE 


ANGULAR MOMENTUM BLOCKADE 


Charge Q = deposit of electric potential energy 


Angular momentum L = deposit of mechanical en- 
ergy 


Electrostatic potential Vraioi'- beyond which there is 
spontaneous discharge 


Angular velocity uj^-^.^^: beyond which the centrifu- 
gal stress results in breaking 


Discharge: decreasing V 


Discharge: decreasing u 


Capacitance C oc area of the plates 


Moment of intertia /in : increasing with spreading 
of the mass distribution about the rotation axis 


Q = CV 


L = /in 


Energy E = CV'^ /2 


Energy E = /in lo"^ /2 


Current i 


Torque r 



TABLE L Mapping of the fundamental quantities and relations relevant to the Coulomb blockade 
to their corresponding quantities for the angular momentum blockade. 

and susceptibility, respectively, which we will now consider. The free energy F is minimized 
by an optimal integer value of / for each value of the applied field B. Hence / forms a 
stepwise function of B with steps at the switching fields B = Bi. The derivative dl/dB 
is a sequence of delta function spikes at the switching fields that should be experimentally 
measurable. We will now consider these phenomena in detail. 
The free energy density can be rewritten as 



F,i\{L,)\=2hl,B) 



F{\{L,)\=2m,B) 
V 



const + aB^ — hBl + cl 



(8) 



forming a piecewise parabolic function, where the following constants were introduced for 

= r/|Aop^,c = ^|Aop^. The critical values 5„ where 
n to I = n + 1 are obtained from the level crossings of the 



convenience: a 



the optimal I switches from / 
functions F{1 = n, B) and F{1 = n + 1, B), which gives 

5„ = ^(2n + l) . 
b 



(9) 



Hence the angular momentum value (L^) = n {2h) is attained for Bn^i < B < Bn- The 
integer n will henceforth denote the integer that minimizes F^, which is given by the integer 
part of (1 + {b/c)B)/2. The magnetization and the susceptibility then become 

dF, 



dB 



-2aB + bn 



(10) 
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and 

X. = ^ = -2a + 6 J] 5(5 - i?„) . (11) 

n 

Except for the non-regular behavior at critical values ([9]) of the field, the magnetic sus- 
ceptibility should have the constant value x = ""i^-^'f^^l^oP- The free energy density 
has cusps joining different parabolic arcs at each critical value for the magnetic field. The 
magnetization has constant jumps inversely proportional to V (connecting regions of linear 
behavior). Finally, the susceptibility has 5-like spikes inversely proportional to V . 



III. BEHAVIOR OF MAGNETIC OBSERVABLES 

Let a; = ^ be the dimensionless area and y = the dimensionless magnetic field. In 
terms of these quantities we have a = -^K^So\Aq\'^x, b = ?7| Aq^-, c = ^| Aq^^, B = j-y = 



2r) X ^ 



— - ~ ■ Then the free energy density, magnetization, and susceptibility are given by 



p - 2^Ef V?/' iVo|Aopyn-n^ 
F^ = const +-K-S,\A,\ ^ (12) 



M^ = -—K—So\Aq\ — y + , (13) 

n 

The magnetization (I13p in terms of the original dimensionful parameters takes the form 



1 I^^^(^IA |2^F , A/'o|Ao|Vb 

^t;(5) = -^K— J Ao — y + —n . (15) 

The second term is responsible for jumps of My at each critical value of the magnetic field. 
All these jumps have the same amplitude and mutual spacing. The first term is due to 
the minimal coupling between the magnetic gauge potential and the order parameter. This 
contains the gradient tensor, and leads to a smooth monotonic variation of the magnetization 
between consecutive jumps. 

The magnitude of the magnetization jumps AMbiockadc associated to the blockade mecha- 
nism increases if the sample area A decreases. In fact, the jump size grows as although 
remaining subleading with respect to the continuous variation AMcont between jumps, even 
for tiny areas A. The first term in Eq. ([T5l) is unaffected by A. 
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FIG. 2. Magnetization AM of the sample (up to an additive constant), plotted in microgauss for 
the example case of area A = (100 nm)^ and thickness d = 100 A, as a function of the magnetic 
field B/Bq, where Bq = Nq{o?c) / {2ri g Ad) ~ 0.7 tesla for such a grain volume. In this plot, the 
field variable B/Bq is centered about a critical level crossing point (i.e. any odd integer) and then 
zoomed by a 10® factor. It is evident how this sample volume entails magnetization jumps ^ l^G . 
The background gradient tensor contribution to the magnetization is by far the leading one, being 
the jumps comparatively micro-scopic (see formula ()16p ): it is represented in units of gauss (with 
no scaling for the axis) in the inset at the top right, for a sample of thickness d = 100 A . 

The ratio between the two contributions to M^i^B) is 



lAMh 



I AM 



cont I 



A ' 



(16) 



where A^.^ = (e^j ~ ^'"^^ A is a critical area independent of any macroscopic 
feature of the SC sample. The magnetization and susceptibility given by Eqs. f lT3|) - f|T^ 
are plotted in Figs. |2}|11 In the plots a cylindrical sample is assumed with dimensions 
V = d ■ A = 100 A X (100 nm)^. The blockade contribution alone to the magnetization is 
depicted in Fig. [3l The dimensionless term K^So\Aq\'^x is proportional to the area A, this 
is the only parameter on which it depends (among the parameters thickness d, area A and 
magnetic field B), and its value affects the plot of the volume susceptibility in Fig. HI 
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FIG. 3. Blockade contribution to the magnetization of the sample, as a function of the magnetic 
field B expressed in tesla, for a cylindrical grain of thickness d = 100 A and area A = (100 nm)^. 
The system displays a staircase of magnetization jumps of the order of 1/U G, which are mutually 
spaced by regular intervals AB ~ 1.4 T. 
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FIG. 4. (Dimensionless) volume susceptibility of the sample, plotted againts the magnetic field B 
for a cylindrical sample of thickness d = lOOA and area A = (100 nm)^; the vertical axis is in units 
of 10^^. There is a background smooth contribution xo to the susceptibility (about —3 • 10~^ for 
this grain) due to the gradient term of the GL functional, which is proportional to the area A of 
the sample, and (5-like spikes above it, due to the angular momentum blockade, which are oc V~^, 
and mutually spaced by regular intervals AB ~ 1.4 T . 
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IV. CONCLUSIONS 



In conclusion we propose the notion of angular momentum blockade in granular d-wave 
superconductors. The effect is due to abrupt changes in the angular momentum of the 
electron liquid in the condensate. The pure d-wave state has Lz = 0. As an out-of-plane 
magnetic field is applied to the sample, a chiral rf'-component is induced. The angular mo- 
mentum of Cooper pairs has to vary accordingly. This process of conversion of momentum, 
from = to a finite value given by the induced d' component, proceeds in a series of 
steps where individual Cooper pairs acquire finite angular momentum. These steps are what 
we call angular momentum blockade. Jumps in the angular momentum of the condensate 
result in jumps of the magnetization, which we calculate. We find that these steps occur 
at well-defined values of magnetic field. We also have defined a "characteristic area" 5*0 (oc 
inverse thickness d~^), providing a scale for the change of induced gaps |Ai| between two 
critical fields causing steps in magnetization, in units of the rf-wave gap amplitude |Ao|. On 
the other hand, we have provided a "critical area" Acr, which is the area scale related to the 
relevant weight between the two distinct contributions entering the magnetization variation 
(within the same magnetic field range considered above): the blockade contribution w.r.t. 
the gradient tensor one. 

As illustrated in Figs. |2|3| a nanoscale grain of typical volume 10^ nm^ should exhibit sig- 
nificant magnetization jumps of the order of 1 microgauss, which are encountered at critical 
magnetic fields whose uniform mutual spacings we estimate being of the order of ~ 1 tesla. 
These results suggest that sensitive magnetization experiments might be able to see these 
jumps. Finally in Fig. H] two different contributions to the magnetic susceptibility are high- 
lighted: the gradient term contribution, which is a constant background value proportional 
to the area A of the sample, and the angular momentum blockade contribution, consisting 
of 5-like spikes oc V^^, and mutually spaced by regular intervals AB (~ 1 tesla for the 
example grain). 
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Appendix A: Relation between the angular momentum and the gaps for the d and 

the d' components 

We refer to the components of the SC wave function, having d^^^yi and d^y symmetry 
respectively, as ipQ and ipi. The angular momentum becomes 

(L,) = {^ij^{e) + ijm\-ihde\MO) + i^i{0)) = 

= -ih J dV '^'^Woi^i cos{2ef - i)oi)l sm{2ef] = -ih J dV - ipoipl] . (Al) 

In terms of respective moduli and phases ipo = \ipo \ exp^iuo), ipi = \ipi \ exp(zz/i), the above 
equation becomes 

(L,) = J t/y[-;i2|^o||V^i|sin(z/o-z/i)] . (A2) 

Because of the global U(l) symmetry we can arbitrarily choose i/q = for this discussion, 
and the free energy is minimized when the relative phase {uq — z^i) = ±7r/2. For a spatially 
uniform wave function, i.e., a spatially uniform order parameter, it is 

(L,) = l^[-;i2|^o||^i| sin(i.o - i^i)] ■ (A3) 
For the minimizing relative phase ±7r/2 this gives 



(L.) = V 



ol 



|2 



|Ai| 



(A4) 



|Ao| 

since |Ai|/|Ao| = |^/'i|/|^/'o|- With a^c the unit cell volume, d the sample thickness, and 
V = Ad the sample volume, the number of Cooper pairs A'sc in the grain is approximately 

V\iJo\' = N,, = g^ , (A5) 

where (7 ~ 0.1 is the high-Tc superconducting fraction. This gives the estimate 

Appendix B: Numerical values 

For allowing better checks, numerical values listed below come with more digits than the 
significative ones. 
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Edges of the SC unit cell 
Coherence length 
Superconducting fraction 
Unperturbed (i^.2_j^2-wave gap 
Density of states 
Fermi energy 
Thickness 

" Characteristic area" 
Modulus of the gradient tensor K 
Dimensionless recurrent combination 
"Critical area" 



a ~ 4A, c ~ lOA 
^ ~ lOA = lO^^cm 
g ~ 10-1 

|Ao| ~ 50 meV = 8.011 x 10"^^ erg 
iVo ~ 2.4 X 10^^ erg"i cm"^ 
Ef ~ 1.98 X 10-12 
c/ = 100 A = 10-^cm 
5-0 = ^ 1.6 X 10-i5cm2 = 16 A' 
~ h-^i'^No ~ 2.156 X 10^3 cm-is-2erg-3 
K^So\Ao\^ ~ 5.68 X 10"^ 
A„ ~ 1.20 X 10-i^cm2 = 1.20 A^ 
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